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At Anderson critical points, the statistics of the two-point transmission Tl for disordered samples 
of linear size L is expected to be multifractal with the following properties [Janssen et al PRB 59, 
15836 (1999)] : (i) the probability to have Tl ~ 1/1/ K behaves as L 9 ^ K \ where the multifractal 
spectrum $(«:) terminates at k — as a consequence of the physical bound Tl < 1; (ii) the 
exponents X(q) that govern the moments T| ~ 1/L x ^ q > become frozen above some threshold: 
X(q > q sa t) = — ${k = 0), i.e. all moments of order q > q sat are governed by the measure of the 
rare samples having a finite transmission (k = 0). In the present paper, we test numerically these 
predictions for the ensemble of L x L power-law random banded matrices, where the random hopping 
Hi,j decays as a power-law (b/\i — j\) a ■ This model is known to present an Anderson transition at 
a — 1 between localized (a > 1) and extended (a < 1) states, with critical properties that depend 
continuously on the parameter b. Our numerical results for the multifractal spectra <!?{,(«) for various 
b are in agreement with the relation <3>(k > 0) = 2 [f(a = d + §) — d] in terms of the singularity 
spectrum f(a) of individual critical eigenfunctions, in particular the typical exponents are related 
via the relation Ktyp(b) = 2(at yp (b) —d). We also discuss the statistics of the two-point transmission 
in the delocalized phase and in the localized phase. 

I. INTRODUCTION 

Since its discovery fifty years ago Anderson localization has remained a very active field of research (see the 
reviews @, H, 0, S S 0> II] ) • One of the most important property of Anderson localization transitions is that critical 
eigenfunctions are described by a multifractal spectrum /(a) defined as follows (for more details see for instance the 
reviews [fj, |8j): in a sample of size L d , the number Ml {a) of points r where the weight If/H^OI 2 scales as L~ a behaves 
as 

M L {a) oc L f ^ (1) 

L — >oo 

The inverse participation ratios (I.P.R.s) can be then rewritten as an integral over a 

Y q (L) = [ d d r\i(j{r)\ 2q ~ / da £ /(q) L- qa ~ L~ T ^ (2) 

J L d J L^oo 

The exponent r(g) can be obtained via a saddle-point calculation in a, and one obtains the Legendre transform 
formula @, 

q =/'(<*) 

T (q) =qa- f(a) (3) 
These scaling behaviors, which concern individual eigenstates ip, can be translated for the local density of states 

p L (E,r) = Y / HE-E n M En (r)\ 2 (4) 

n 

as follows : for large L, when the L d energy levels become dense, the sum of Eq. [5] scales as 

p L (E,r) oc L d \yj E (r)\ 2 (5) 
and its moments involve the exponents r(q) introduced in Eq. [5] 



[p L (E,r)]« with A(q)=T(q)-d(q-l) (6) 

These notions concern one-point functions, and it is natural to consider also the statistics of two-point functions. 
In particular, a very interesting observable to characterize Anderson transitions is the two-point transmission Tl 
when the disordered sample of size L d is attached to one incoming wire and one outcoming wire [3, l2a | : (i) it 
remains finite in the thermodynamic limit only in the delocalized phase, so that it represents an appropriate order 
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parameter for the conducting/non-conducting transition; (ii) exactly at criticality, it displays multifractal properties 
in direct correspondence with the multifractality of critical eigenstates, i.e. it displays strong fluctuations that are 
not captured by more global definitions of conductance. More precisely, as first discussed in [9[ for the special case 
of the two dimensional quantum Hall transition, the critical probability distribution of the two-point transmission Tl 
takes the form 

Prob (T L ~ L~ K ) dT cx L^dn (7) 

L—*oo 

and its moments involve non-trivial exponents X(q) 

T£~ [ dnL^- qK cx L~ X W (8) 

As stressed in Q, the physical bound Tl < 1 on the transmission implies that the multifractal spectrum exists only 
for k > 0, and this termination at k — leads to a complete freezing of the moments exponents 

X(q)=X(q aat ) for q > q sat (9) 

at the value q sa t where the saddle-point of the integral of Eq. [H] vanishes n(q> q sa t) = 0. It is very natural to expect 
some relation between the two multifractal spectra f(a) and and the possibility proposed in |9( is that before 

the freezing of Eq. [5] occurs, the transmission should scale as the product of two independent local densities of states 
(Eq. 

X(q) = 2A(q) for q < q sat (10) 

We refer to @ for physical arguments in favor of this relation. Equations [9] and [10] for the moments exponents are 
equivalent to following relation between the two multifractal spectra 

*(/s>0)=2[/(a = d+|)-d] (11) 

In this paper, our aim is to test numerically these predictions for the statistics of the two-point transmission Tl 
at the critical points of the Power-law Random Banded Matrix (PRBM) model, where one parameter allows to 
interpolate continuously between weak multifractality and strong multifractality. We will also discuss the statistics 
of the two-point transmission off criticality. 

The paper is organized as follows. In section [HI we introduce the PRBM model and the scattering geometry used 
to define the two-point transmission. In Section IIII1 we present our numerical results concerning the multifractal 
statistics of the two-point transmission at criticality. We then discuss the statistics of the two-point transmission in 
the localized phase fSection lTV)) and in the delocalized phase (Section|V| respectively. Our conclusions are summarized 
in section IVT1 The appendices A and B contain more details on the numerical computations. 



II. MODEL AND OBSERVABLES 



Beside the usual short-range Anderson tight-binding model in finite dimension d, other models displaying Anderson 
localization have been studied, in particular the Power-law Random Banded Matrix (PRBM) model, which can be 
viewed as a one-dimensional model with long-ranged random hopping decaying as a power-law (b/r) a of the distance 
r with exponent a and parameter b (see below for a more precise definition of the model). The Anderson transition 
at a = 1 between localized (a > 1) and extended (a < 1) states has been characterized in [l(| via a mapping onto 
a non-linear sigma-model. The prop erties of the critical points at a = 1 have been then much studied, in particular 
the statistics of eigenvalu es fill [l2[ [r| , and the multifractality of eigenfunctions [3, [ID, [TH, [l?], [H, 13 , including 
boundary multifractality [20[. The statistics of scattering phases, Wigner delay times and resonance widths in the 
presence of one external wire have been discussed in [2l|, [2J| ■ Related studies concern dynamical aspects [23[ , the case 
with no on-site energies (24j, and the case of power-law hopping terms in dimension d > 1 [25|, [2(| [2tJ • In this paper, 
we consider the PRBM in a ring geometry (dimension d = 1 with periodic boundary conditions) in the presence of 
two external wires to measure the transmission properties. 



3 



A. Power-law random banded matrices with periodic boundary conditions 

We consider L sites i = 1,2,..L in a ring geometry with periodic boundary conditions, where the appropriate 
distance r^j between the sites i and j is defined as [3| 

r u=~ iim [ t ( 12 ) 



**» 7T \ L 

The ensemble of power-law random banded matrices of size L x L is then defined as follows : the matrix elements 
Hij are independent Gaussian variables of zero-mean Hij = and of variance 



H h = (13) 

The most important properties of this model are the following. The value of the exponent a determines the localization 
properties [lOj : for a > 1 states are localized with i nteg rable power-law tails, whereas for a < 1 states are delocalized. 
At criticality a = 1, states become multifractal flil Il5l Aldi and exponents depend continuously of the parameter 
b, which plays a role analog to the dimension d in short-range Anderson transitions : the limit b 3> 1 corresponds 
to weak multifractality ( analogous to the case d = 2 + e) and can be studied via the mapping onto a non-linear 
sigma-model [10|j , whereas the case b <C 1 corresponds to strong multifractality ( analogous to the case of high 
dimension d) and can be studied via Levitov renormalization [lil . [3J| . Other values of b have been studied numerically 



B. Scattering geometry used to define to two-point transmission 




FIG. 1: (Color on line) The ensemble of power-law random banded matrices of size L x L can be represented as a ring of L 
sites, where the matrix element Hij between the sites i and j is a Gaussian variable of zero-mean Hij = and of variance 
given by Eq. [13] in terms of the distance r^j of Eq. [T2j In this paper, we study the statistics of the Landauer transmission 
T L = \t\ 2 when an incoming wire is attached at the site L/2 and an outgoing wire is attached at the site L (see text for more 
details) . 

In quantum coherent problems, the most appropriate characterization of transport properties consists in defining a 
scattering problem where the disordered sample is linked to incoming wires and outgoing wires, and in studying the 
reflection and transmission coefficients. This scattering theory definition of transport, first introduced by Landauer 
[29l ] , has been often used for one-dimensional systems [3(J [32j and has been generalized to higher dimensionalities 
and multi-probe measurements (see the review [HI). In the present paper, we focus on the Landauer transmission 
for the scattering problem shown on Fig. [1] an incoming wire is attached at the site L/2 and an outgoing wire is 
attached at the site L. We are thus interested into the eigenstate \ip > that satisfies the Schrodinger equation 



H\i/> >= E\tp > 



(14) 
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inside the disorder sample characterized by the random Hij, and in the perfect wires characterized by no on-site 
energy and by hopping unity between nearest neighbors. Within these perfect wires, one requires the plane-wave 
forms 

*Pout(x>x L ) =t e lk ^- x ^ (15) 

These boundary conditions define the reflection amplitude r of the incoming wire and the transmission amplitude t 
of the outgoing wire. The Landauer transmission 

T^| 2 = i-kf (i6) 

is then a number in the interval [0, 1]. More details on the numerical computation of the transmission in a given 
sample are given in Appendix |A"1 

To satisfy the Schrodinger Equation of Eq. [TJ] within the wires with the forms of Eq. [TS1 one has the following 
relation between the energy E and the wave vector k 

E = 2cosk (17) 

To simplify the discussion, we will focus in this paper on the case of zero-energy E — (wave- vector k = 7r/2) that 
corresponds to the center of the band. 

In the following, we study numerically the statistical properties of the Landauer transmission T for rings of size 
50 < L < 1800 with corresponding statistics of 10. 10 8 > n s (L) > 2400 independent samples. For typical values, the 
number n s (L) of samples is sufficient even for the bigger sizes, whereas for the measure of multifractal spectrum, we 
have used only the smaller sizes where the statistics of samples was sufficient to measure correctly the rare events. 



III. STATISTICS OF THE TWO-POINT TRANSMISSION AT CRITIC ALITY (a = 1) 

As recalled in the Introduction, the two-point transmission is expected to display multifractal statistics at 
criticality @. We first focus on the scaling of the typical transmission before we turn to the multifractal spectrum 
and the moments of arbitrary order. 



A. Typical transmission at criticality (a — 1) 

As discussed in [t| [28[ , the typical transmission 



L 

is expected to decay at criticality with some power-law 



typ _ JnT L 



1 



oc — — (19) 

where the exponent K typ is directly related via the relation 

K t y P = 2{at yp - d) (20) 
to the typical exponent ctt yp that characterizes the typical weight of eigenfunctions 

(21) 

This typical value at yp corresponds to the maximum value f(at yp ) = d of the multifractal spectrum f(a) introduced 
in Eq. Q] . (Note that in 0, [2^|), n typ is denoted by X t and a typ by ctQ. Here we have chosen to use the explicit 
notation 'typ' for clarity). 

For the PRBM considered here, the dimension is d = 1, and critical exponents depend continuously on b. We show 
on Fig. [2] (a) the hiT L yp as a function of InL : the slopes allows to measure the exponents K typ {b) 



\nT L vp {a^l,b) = \nT L {a=l,b) oc -n typ {b)\nL (22) 

L — >OG 
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FIG. 2: (Color on line) Scaling of the typical transmission T l ^ v at criticality a = 1 for various values of the parameter b : (a) 
]nT^" p = InTi, as a function of InL for various values of b = 4 (V), 1 (O),0.5 (□),0.25 (0),0.1 (a),0.05 (>) (the two other 
values b = 10 and b — 0.01 we have studied are not shown here for clarity) : the slope yields the typical exponent Ktyp(b) of 
Eq. 1221 (b) critical exponent Kt y p{b) as a function of b. 



b 


b^O 


0.01 


0.05 


0.1 


0.25 


0.5 


1 


4 


10 


b — » +oo 


Ktyp(b) 


2 


1.92 


1.64 


1.33 


0.77 


0.46 


0.25 


0.08 


0.04 





a typ (b) = 1 + SMEW 


2 


1.96 


1.82 


1.66 


1.38 


1.23 


1.12 


1.04 


1.02 


1 



TABLE I: Critical exponents as a function of b : (i) the exponent «t ap (6) characterizes the typical transmission at criticality 
(see Eq. 122ft . (ii) the corresponding value of the typical exponent a typ = 1 + Kt yP (b)/2 (see Eq. I20[) for the weight ip 2 {r) of 
eigenfunctions (see Eq. I21|l 



On Fig. [2] (b), we show how the exponent Kt yp {b) depends on b. The values Kt yp (b) we have measured are listed in 
Table HI together with the corresponding values of ott yp (b) obtained via Eq. [50] : these values of at yp (b) are compatible 
with the values of the maxima of the multifractal spectrum f(a) of critical eigenstates measured in [l4[ (see Fig. 2 
and Fig. 6 of ;i4j) and in [H (see Fig. 2 and Fig. 3 of [H). 

The two limits 6^1 and b <c 1 can be understood as follows. The case 6^1 corresponds to very weak 
multifractality with the typical exponent cttyp ~~ > 1 [3 El ■ Equation [50] yields that the critical exponent n typ of the 
typical transmission becomes arbitrary small in the limit b — > +oo 

+°°) -* ( 23 ) 

The opposite limit 6^1 corresponds to very strong multifractality with the typical exponent at yp —* 2 
Equation [2^1 thus yields 

iHy P (b -> 0) -> 2 (24) 



B. Multifractal spectrum $&(«;) with termination at k = 

As recalled in the introduction, the statistics of the two-point transmission is expected to be multifractal at crit- 
icality, as a consequence of the multifractal character of critical eigenfunctions @. For the PRBM, the multifractal 
spectrum $(,(«) of Eq. [7] will depend continuously on the parameter b 

Prob (T L - L~ K ) dT oc L^ck (25) 

Since it describes a probability, the multifractal spectrum satisfies $;,(«;) < 0, and the maximal value $&(/«) = is 
reached only for the typical value k — K typ (b) discussed above 



®b{Kt VP (b)) = 



(26) 
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FIG. 3: (Color on line) Multifractal statistics of the critical two-point transmission for three values b — 1, b — 0.1 and 
b = 0.01 : (a) the multifractal spectra ^b(n) reach their maximum $b(K typ (6)) = at the typical values Kt y p(b) and exactly 
terminate at finite values "^(k = 0). (b) the corresponding moments exponents X(,((j>) become completely frozen for q > q sa t '■ 
X b (q > qsat) = = 0). 

The relation of Eq. 1201 between the two typical exponents K tvp (b) and at yp (b) is expected to come from the more 
general relation of Eq. ITT1 between the two multifractal spectra and fb{a) 

$ b (K> 0) = 2 [f b (a = 1 + |) - l] (27) 

An essential property of the spectrum ${,(k) is that it exists only for k > as a consequence of the physical bound 
Tl < 1, so that it terminates at k = at the finite value 

$ h ( K = 0) = 2 [f b (a = 1) - 1] (28) 

To measure numerically the multifractal spectrum <£>;,(«;), we have used the standard method based on g-measures 
of Ref. [35| (see more details in Appendix [B]) . To show how the parameter b allows to interpolate between weak 
multifractality and strong multifractality, we compare on Fig. [3] (a) the multifractal spectra $;,(k) for the three values 
b = 1, b = 0.1 and b = 0.01. For instance for the value b — 0.1, the termination value we measure $fc(K = 0) ~ —0.58 
is in agreement via Eq. [55]with the value fb{a — 1) ~ 0.71 of Fig. 6 of [3| and Fig. 2 of (l9j . 

C. Freezing transition of the moments exponents Xb(q) 
As usual, the multifractal statistics of Eq. [21] has for consequence that the moments of arbitrary order q 

Tf ~ / dnL* b ^- qK oc L~ Xb ^ (29) 

J L^oo 

are governed by non-trivial exponents X b {q) that can be obtained via the saddle-point calculation 

- X b (q) = max [$ 6 ( K ) - nq] (30) 

K>0 

As long as the saddle-point value satisfies n{q) > 0, Xb(q) can be obtained via the usual Legendre transform formula 

q 

X b (q) =Kq-$ b (K) (31) 
However above some threshold q sa t, the saddle-point value will saturate to the boundary value 

«(g > q sat ) = o (32) 
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and the exponent X{q) will saturate to the value 

X b (q > q sat ) = X h (q sat ) = -$ 6 (k = 0) = 2(1 - f b (a = 1)) (33) 

This freezing phenomenon of X{q) at q sa t predicted in [9( means that all moments of order q > q sat are dominated 
by the rare events corresponding to a finite transmission T ~ 1, whose measure behaves as L*( K=n ). 

We show on Fig. [3] (a) the moments exponents Xt(q) for the three values 6 = 1,6 = 0.1 and 6 = 0.01. For instance 
for b = 0.1, the freezing value X b {q > q S at) ~ 0.58 corresponds to the termination value = 0) ~ —0.58 of Fig. [3] 

It turns out that for Anderson transitions, a special symmetry of the multifractal spectrum /(a) has been proposed 
(see [Til |36[ and references therein) that relates the regions a < d and a > d via the relation f(2d — a) = f(a) + d — a. 
This symmetry then fixes the value of q sat where n{q sa t) = or equivalently a(q sat ) = d to be exactly 

Qsat = \ (34) 

Numerically, it is difficult to measure precisely the value q sa t where the exponents X (q) become completely frozen 
as a consequence of finite-size corrections around this phase transition point for the X(q), as already found for the 
quantum Hall transition in Q. However Fig. [3] (b) shows that in the limit of strong multifractality (b = 0.01), the 
numerical saturation value is not far from the theoretical prediction of Eq. 1341 

IV. STATISTICS OF THE TWO-POINT TRANSMISSION IN THE LOCALIZED PHASE (a > 1) 

A. Typical transmission in the localized phase 



-5 




-35 



In L 



FIG. 4: (Color on line) Typical two-point transmission in the localized phase a > 1. The power-law decay of Eq. [35]is checked 
here for the case b = 0.1 and the two values a = 1.4, a = 1.8 : the slopes of this log-log plot are of order 2a ~ 2.8 and 2a ~ 3.6. 
For comparison, we also show the critical data for a — 1 of slope n tyP (b = 0.1) ~ 1.33. 

In usual short-range models, the localized phase is characterized by exponentially localized wavefunctions, whereas 
in the presence of power-law hoppings, localized wavcfunction can only decay with power-law integrable tails. For 
the PRBM, it is moreover expected that the asymptotic decay is actually given exactly by the power-law of Eq. [T^] 
for the hopping term defining the model [l(| : IV'Mltyp ~ l/r 2a . As a consequence in the localized phase a > 1, one 
expects the typical decay 



rptyp 



(35) 



As shown on Fig. [4j we have checked this power-law decay of the typical transmission for the case b = 0.1 and the 
two values a — 1.4, a — 1.8. 
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In T InT 

FIG. 5: (Color on line) Histogram of of the two-point transmission in the localized phase a > 1 : example for b = 0.1 and 
a = 1.4 : (a) histograms Pl of In Tz, for the four sizes L — 100,200,400,800 (the histograms for bigger sizes are more noisy 
and are not shown for clarity) (b) same data in log scale to exhibit the tails discussed in the text. 



B. Histogram of (In Tl) in the localized phase (a > 1) 

We show on Fig. [5]the histograms Pl of InTf, for the four sizes L = 100, 200, 400, 800 in the localized phase a = 1.4 
for the case b = 0.1 : these data seem to indicate that as L grows, the probability distribution shifts to the left while 
keeping a fixed shape. This means that the relative variable u = (lnT^ — \nT L yp ) with respect to the typical value 
discussed above (see Eq. 135]) remains a finite variable as L — > +oo. In addition, the left tail is governed by the 
exponent a ~ 0.5 

hxP^(u=(lnT L -lnTly p )) ~ \u (36) 
or equivalently after the change of variable T — = e u 

pl ^ T = ^K*oT^ (37) 

These properties can be understood by the following simple argument. In the localized phase a > 1, one may 
assume that for large L, the transmission Tl is dominated by the direct hopping term Hl/2,l (see Fig. 1) 

2 

T L (a >1) oc |Jf w | 2 « ^ (38) 

where x is a Gaussian variable of zero mean and variance unity (see Eq. I13p . In particular, its probability density is 
finite near the origin P(x = 0) > 0. The change of variable 



1 L 



then yields the power-law of Eq. [371 



V. STATISTICS OF THE TWO-POINT TRANSMISSION IN THE DELOCALIZED PHASE (a < 1) 

A. Typical transmission in the delocalized phase 

In the delocalized phase, the eigenfunctions are not multifractal anymore, but monofractal with the single value 
&deioc = d for the weight As a consequence, the typical transmission is expected to remain finite as L — > +oo Q (in 
Eq. [201 the case a typ = d yields K typ = 0) 

T L yp cx T^ p >0 (40) 

L— *OQ 

The two-point transmission is thus a good order parameter of the transport properties Q . 
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B. Histogram of (In Tl) in the delocalized phase a < 1 




FIG. 6: (Color on line) Histogram of of the two-point transmission in the delocalized phase a < 1 : example for b = 1 and 
a — 0.75 (a) the histograms of IuTl for the three sizes L = 100,200,400 coincide (the histograms for bigger sizes are more 
noisy and are not shown for clarity) (b) same data in log scale to exhibit the tail (see Eq. 141 p 



We find that the whole probability distribution Pt(lnT) converges for large L towards a fixed distribution 
•Pjxf (haT). As an example, we show on Fig. [H] the histograms of (InTj,) for three sizes L = 100,200,400 con- 
cerning the case 6=1 and a = 0.75. These histograms stops at InT = as a consequence of the bound T < 1. In the 
region of very small transmission InT — > — oo the log-log plot of Fig. [5] (b) indicates the same form as in Eq. [351 



or equivalently after a change of variable 



lnP^ oc (lnT) cx J InT (41) 

In T — > — oo z 



t el ° c (T) T « Q ^ (42) 



As explained around Eq. 1391 this power-law behavior simply means that the transmission can be written as the square 
T ~ x 2 of some random variable x that has a finite probability density at the origin P(x = 0) > 0. This means that 
all negative moments of order q < —1/2 actually diverge in the delocalized phase 



l 

deloc 



dTT q P^ eloc (T) ~ +oo (43) 

g<-l/2 



VI. CONCLUSIONS AND PERSPECTIVES 



In this paper, we have studied numerically the statistics of the two-point transmission Tl as a function of the size 
L in the PRBM model that depends on two parameters (a, b) : 

(i) in the delocalized phase (a < 1), we have found that the probability distribution of Tl converges for L — > +oo 
towards a law P^ loc (T) presenting the power-law of Eq. [42] for T — > 0. 

(ii) in the localized phase (a > 1), we have found that the probability distribution P l oc (\uTl) keeps a fixed shape 
around the typical value lnT* yp = lnT^ as L grows, and the typical value decays only as the power-law of Eq. 
[35l as a consequence of the presence of power-law hoppings. 

(hi) exactly at criticality (a = 1), the statistics of the two-point transmission Tl is multifractal : we have measured 
the multifractal spectra $&(«;) as well as the moments exponents Xb(q) for various values of the parameter b that 
allows to interpolate between weak multifractality and strong multifractality. We have tested in detail the various 
expectations of Ref. Q concerning the termination of $&(«;) at k = 0, the freezing of Xb(q) above some value q > q sa t, 
and the relations with the singularity spectrum fb(ot) of individual critical eigenstates. 

To finish, we should stress that the relation of Eq. [TT] relates the transmission between two 'bulk points' to the 
'bulk multifractal spectrum' f{a). In localization models, it is however more usual to attach leads to the boundaries 



10 



of the disordered sample : then the statistics of the two-point transmission is related to the 'surface multifractal 
spectrum' as will be discussed in more details elsewhere (37]] • We will also discuss in 37] the statistical properties of 
the transmission for various scattering geometries involving a large number of wires. 

APPENDIX A: COMPUTATION OF THE TWO-POINT TRANSMISSION IN EACH SAMPLE 

To compute the transmission of a given sample via Eq. I16i we have to solve the Schodinger problem of Eq. 1 141 with 
the scattering boundary conditions of Eq. [15] This can be decomposed in two steps as follows. 

1. Recursive elimination of the 'interior sites' 

The first step consists in the iterative elimination of the 'interior sites', i.e. of all the sites not connected to the 
external wires (see Fig. Q]). To eliminate the site io, one uses the Schodinger Eq. [14] projected on this site 

£ty(*o) = H io , io ^{i Q ) + ]T H i0 ji>(j) (Al) 

3 

to make the substitution 

*Kio) = w \ E H io,M3) (A2) 

in all other remaining equations. Then from the point of view of remaining sites, the hoppings are renormalized 
according to 

This procedure is stopped when the only remaining sites are the two sites L/2 and L connected to the external wires 
(see Fig. [1} : the three real remaining parameters are the renormalized hopping H L / 2 L an< i th e two renormalized 
on-site energies H l/2 ,l/2, 

2. Effective scattering problems for the two boundary sites 

The second step consists in solving the scattering problem for the two boundary sites L/2 and L connected to the 
external wires (see Fig. [T]) with the renormalized parameters obtained above. The Schodinger Eq. [14] projected on 
the boundary sites L/2 and L read 

Eip(x L/2 ) = H l/2 ^ l/2 i/j(x l/2 ) + ip(x L/2 - 1) + H L/2tL ip(x L ) 
Ei>{x L ) = H L .Lip{x L ) + ip{x L + l) +H LtL/2 il>(x L/a ) (A4) 
The boundary conditions of Eq. [15] fixes the following ratio on the outgoing wire 

+ = e , k (A5) 
lp{XL) 

The following ratio 

concerning the incoming wire can be then computed in terms of the three real renormalized parameters 

E - ( h l . l + e k 



H 2 

R = E- H L/2 , L/2 ,_ L ' L/2 - (A7) 



The reflexion coefficient r of Eq. [15] is then obtained as 

R - e~ ik 



e ik - R 

yielding the transmission of Eq. 1161 



(A8) 
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APPENDIX B: COMPUTATION OF THE MULTIFRACTAL SPECTRUM OVER THE SAMPLES 

To measure numerically the multifractal spectrum <I>(k) of Eq. [7]that characterizes the statistics of the transmission 
T L over the samples of size L, we have used the standard method based on q- measures of Ref. (35[. More precisely, 
for various sizes L, we have measured the transmission T^(i) for a number n s {L) of independent samples (i). Then 
for various values of q, we have computed the moments of Eq. [8] 

^^gW (Bl) 

to extract the moments exponents X(q) as the slope of the log- log plot 

ln(rf) oc -X(q)lnL (B2) 
We have also computed the auxiliary observables 

n B (L) 

E [T L {i)] q (-^T L (i)) 

E p^)]' 

and 

F L (q) = q K L (q)+\n(jfj (B4) 

to obtain n(q) and as the slopes of 



and 



K L (q) oc n(q)\nL (B5) 

L — >co 



F L (q) oc ${ K (q))]nL (B6) 

L — >oo 



This yields a parametric plot in q of the multifractal spectrum $(k) : on Figured] (a), each circle of coordinate 
(k, $(«;)) corresponds to a value of q. 
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